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Abstract 

We investigate exact eigenstates of tight-binding models on the planar rhom- 
bic Penrose tiling. We consider a vertex model with hopping along the edges 
and the diagonals of the rhombi. For the wave functions, we employ an ansatz, 
first introduced by Sutherland, which is based on the arrow decoration that 
encodes the matching rules of the tiling. Exact eigenstates are constructed for 
particular values of the hopping parameters and the eigenenergy. By a gener- 
alized ansatz that exploits the inflation symmetry of the tiling, we show that 
the corresponding eigenenergies are infinitely degenerate. Generalizations and 
applications to other systems are outlined. 

PACS numbers: 71.23.Ft, 05.60.+W, 71.23.An, 71.30.+h 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The discovery of quasicrystals by Shechtman et alB stimulated wide interest in the physics 
of these materials which are intermediate between periodic and random structures. Besides 
icosahedral quasicrystals, which are aperiodic in all three dimensions of space, also periodi- 
cally layered structures with planar quasiperiodic order and non-crystallographic rotational 
symmetries were found, comprising dodecagonal,i decagonal,! and octagonal^ phases with 
twelve-, ten-, and eightfold symmetry, respectively. Although the fundamental question 
"where are the atoms?", raised e.g. in Ref. |, has only been answered partially to date, 
most structure models of quasicrystals are based on two- or three-dimensional quasiperiodic 
tilings or their disordered versions (random tilings). 

An important and exciting problem in condensed-matter physics is whether the quasiperi- 
odic structure leads to new and unexpected physical properties. In particular transport 
properties, as for instance electric or heat conductance, are strongly effected by the non- 
periodic order. Indeed, many quasicrystalline alloys are characterized by very high values of 
electric resistivity, by a negative temperature coefficient of resistivity, and by a low electronic 
contribution to the specific heat which points to a small density of states at the Fermi energy. 
It is difficult to explain these striking features, because a rigorous theory of the electronic 
structure of quasiperiodic materials does not exist. For want of a simple analogy of Bloch 
theory for quasicrystals, one either carries out numerical calculations for as large clusters 
as possible, or one tries to make exact statements about the electronic wave functions in 
simple models. 

We consider tight-binding models on the two-dimensional rhombic Penrose tiling.! Our 
models are so-called vertex models because we locate the atoms at the vertices of the tiling. 
Interactions are taken into account only between neighboring vertices connected by edges 
or by diagonals of the rhombi. In our calculations, we restrict ourselves to a single s-type 
atomic orbital per vertex. This makes the transfer integrals tij independent of the angular 
orientation and leads to the following Hamiltonian 



where \i) denotes a Wannier state localized at vertex i, and Ei are on-site energies. For 
the hopping integrals t^, we choose five different values 1, di, d2, d^, d^, depending on the 
distance of the vertices i and j, see Fig. |l]. Here, tij = 1 for vertices connected by an edge of 
the tiling, tij = di (^2) for the long (short) diagonal of the 'fat' rhombus, and tij = d^ (^4) 
for the long (short) diagonal of the 'thin' rhombus, respectively. 

As the Penrose tiling is arguably the most popular among the quasiperiodic tilings, 
it is not surprising that tight-binding models defined on the Penrose tiling have been in- 
vestigated rather thoroughly. Besides the vertex model,0"Ei the so-called center model was 
consideredjn^^ where atoms are located in the center of the rhombi, and hopping may occur 
between adjacent tiles — this is nothing but a vertex model on the dual graph of the Penrose 
tiling. However, most results rely on numerical approaches, and only few exact results on 
the spectrum of the tight-binding Hamiltonian are known. In particular, so-called 'confined 
states' have been investigated in detail, both for the vertexifl and the center model.il These 
are infinitely degenerate, strictly localized eigenstates corresponding to a particular value of 
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the energy, which occur as a consequence of the local topology of the tihng, see also Ref. [T8 . 
Furthermore, for a Hamiltonian (|1 . 1| ) with particular on-site energies Ei chosen according to 
the vertex type at site i, the exact self-similar ground state could be constructed.^ Based 
on the same idea, several non-normalizable eigenstates of the center model and their multi- 
fractal properties were obtained exactly.^ These solutions, restricted to special values of the 
hopping integrals, were derived from a suitable ansatz for the eigenfunctions. According to 
this ansatz, the wave function at a site depends only on its neighborhood and on a certain 
integer number associated to the site, a 'potential', which is derived from the matching rules 
of the Penrose tiling.0 

In this paper, we apply the same ansatz to the vertex model on the Penrose tiling. The 
solution is more complicated than for the center model, where the coordination number 
(i.e., the number of neighbors) is always equal to four, whereas for the vertex model it 
varies between three and seven, or between six and fourteen if we include neighbors along 
diagonals, respectively. For suitably chosen transfer integrals, we derive exact eigenstates of 
the Hamiltonian ( |1 . 1|) and analyze their multifractal behavior. As observed for the center 
model,ii we find that these states are infinitely degenerate, i.e., for fixed value of the en- 
ergy the eigenfunctions still involve one free parameter. In order to show this, we need to 
generalize the ansatz exploiting the inflation symmetry of the Penrose tiling. 

Our paper is organized as follows. In the subsequent section, we discuss the labeling of 
the rhombi with two kind of arrows and the associated potentials. In Sec. |Tl|, we introduce 
the ansatz for the wave function and solve the tight-binding equations for two cases, the first 
one with Ei = and di ^ 0, and the second one with di = but various on-site energies. A 
generalized ansatz, based on the inflation symmetry of the tiling, is considered in Sec. |IV|. In 
Sec. 0, we perform a fractal analysis of the wave functions, i.e., we calculate the generalized 
dimensions. Finally, we conclude in Sec. |V1[ 



II. EDGE-LABELING AND POTENTIALS 

Following de Bruijn,@ we mark the rhombi with single and double arrows as shown in 
Fig. 1^. The matching rules require that arrows on adjacent edges match. Fixing a certain 
site O as the origin, we assign to a site i two integers n{i) and m{i) which count the number 
of single and double arrows, respectively, along an arbitrary path connecting the origin O 
and site i. This is well-defined because, along any closed path, the total number of single 
and double arrows vanishes, as can be seen from Fig. ^. We refer to these integers n{i) and 
■m(i) as 'potentials' at site i because they are integrals of the two vector fields defined by 
the arrows. The distributions of the potentials are rather irregular and show the following 
properties: 

• The single-arrow potential n{i) is directly related to the sum t{i) G {1,2,3,4} of the 
five-dimensional indices denoting the translation class of the site i. It takes only two 
values: n{i) = if t{i) G {2,3} and n{i) = 1 if t{i) G {1,4} (provided the origin has 
translation class t(0) G {2,3}). 

• The double-arrow potential m{i) is unbounded. Its distribution on a finite patch is 
shown in Fig. |^. For a detailed discussion of the distribution, see Ref. |Tp . 
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The potential m{i) is the key ingredient in the construction of exact eigenf unctions of tight- 
binding Hamiltonians on the Penrose tihng. 

III. SOLUTIONS OF THE TIGHT-BINDING EQUATIONS 

We want to construct solutions of the tight-binding equations 

Y.tij(j)j = iE-ei)<Pi, (3.1) 
j 

where we sum over all neighbors j of the site i. As an ansatz, we demand that the wave 
function amplitude 0j at site i depends solely on the vertex type z/(i) G {1, 2, . . . , 8} and on 
the potential m{i). This leads to the following ansatz 

0, = . (3.2) 

The eight vertex types of the Penrose tiling are shown in the top row of Fig. ^. The 
corresponding eight amplitudes A^, and /5 are parameters. 

A. The case = 

For simplicity, we first concentrate on the case with on-site energies Ei = 0. With the 
ansatz ( p.2| ), the infinite set of equations ( |3.1| ) reduces to a finite set comprising as many 
equations as there are second-order vertex types in the tiling. By a second-order vertex 
type we mean the neighborhood of a site up to its second coordination zone. There are 31 
different second-order vertex types in the Penrose tiling, these are shown in Fig. ^, grouped 
together according to the first-order vertex type of the central site given in the top row. 
Thus, we have 31 linear equations in the 14 variables A^, = 1, . . . ,8), j3, di, d2, d^, ^4, 
and E. As it is straightforward to derive the equations from the second-order vertex types 
of Fig. ^ we refrain from listing them here. Instead, we consider as an example only the 
second-order vertex types in the first column of Fig. ^, which we show again in Fig. ^ (rotated 
by 90 degrees) together with the corresponding values of the potential m{i). This yields the 
following four equations 

EA, = d^Asp + 2{A2 + A3) + di{A2 + + A7)/3-^ 
= da^P + 2(^2 + A3) + di{A2 + A, + A,)p-' 
= d^Asp + 2{A2 + A3) + di{A2 + A-r + Ag)/?"^ 
= d^Asp + 2{A2 + A3) + di{A2 + + A7)p-^ 

(3.3) 

two of which (the first and the third) are identical because the corresponding patterns are 
mirror images of each other. 

At first sight, as the number of variables, 14, is much smaller than the number of equa- 
tions, 31, one might expect that only the trivial solution (0j = 0) exists. However, this is 
not the case, for suitably chosen values of the hopping parameters di, d2, d^, d^, and the 
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energy E, non-trivial solutions exist, because the equations are not independent. To see 
this, note that the second-order vertex types within one column of Fig. ^ differ only slightly 
from each other, which means that the corresponding equations are also very similar as can 
be seen in the example (|3.3| ). Thus, they can be substantially simplified by subtraction. For 
example, the differences between the equations in (|0|) result in the single equation 



^1(^5-^7) = (3.4) 

which implies A<^ = A-j (unless d\ vanishes). From the analogous equations for the other 
vertex types, it turns out that the amplitudes Ay{i-) depend only on the translation class t{i) 
of the site rather than on its specific vertex type This means 

A^ = A^ = A^ (te{l,4}) 

A2 = A^ = A^ = A, = As it G {2,3}) . (3.5) 

With this, all equations corresponding to second-order vertex types with the same central 
vertex reduce to a single equation, and one is left with the following eight equations 

EAx = diA2f3 + AA2 + 3^1^2/3"^ 

= 5^2 + MiA^lS-^ 

= 2d4A2(3 + 3^2 + diA^r^ 
EA2 = {diA, + 2A2)I3 + 2Ai + 2(^2 + 4)^2 + ^2/?"' 

= Ai + 2d2A2 + {2A2 + d^A^)p-^ 

= {3diAi + 5A2)f3 + 2A, + 44^2 

= {AdiAi + 5A2)(3 + Ai + 2d^A2 

= (5diAi + 5^2)/? (3.6) 

for central vertices of type 1, 4, 6, and 2, 3, 5, 7, 8, respectively. 

For this system of equations, we obtain three sets of non-trivial solutions, expressed in 
terms of the parameter /3 which may be chosen freely. Here, we introduce 
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[3 ±(3-^ 



(3.7) 



to abbreviate the formulae below. 

Solution (1): The wave function has the form 



(l-2/52)/5™W fort(2)G{l,4} 

l^m{i)+i for t(z) G {2,3} ^ ' 



for transfer integrals and energy given by 

di = \h_ d2 = -jb^ + b 

is = + \ 

5 
2 

Solution (2): Here, 



2"- ^2. 
ds = -jb^ + \hZ^ di = b- 

E = (3.9) 
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^m(i)+i for t(i) G {1,4} 
f^m(i) fort(i)G{2,3} 



(3.10) 



with 



E = 56+ (3.11) 



Solution (3): Finally, 



for t(«) G {2,3} ^ ^ ^ 



where 



c?i = -(| + /5')&+ d2 = \b+-\b-^ 

4 = -(i + r')&+ rf4 = -&+ 

E = §6+ (3.13) 

For each of these solutions, there exists an additional solution for a slightly generalized 

ansatz 

0^ = i.W,tW/3"(*) (3.14) 

that involves the translation class t{i) at site i. Note that for each vertex type there are only 
two possible values, t{i) e {1, 4} for iy{i) e {1, 4, 6} and t{{} e {2, 3} for z/(i) e {2, 3, 5, 7, 8}, 
which were not distinguished in our previous ansatz, see Eq. (|3.5| ). The wave functions differ 
from the solutions given in Eqs. ( |3.8| ), ( |3.10| ), and ( p.l2| ) only by an alternating sign which 



depends on the translation class 

0, = (-l)*«0i (3.15) 

and by a sign change in the parameters, i.e., di = —di, d2 = —d^, d^ = —d^, d^ = —d^, 
and E = —E. Note that the two models differing by this sign change are not trivially 
related, because the hopping parameter along the edges of the tiles does not change its sign 
— it is always equal to 1. These six solutions exhaust all non-trivial solutions in terms of 
the ansatz ( |3.14| ), but we note that for a given set of hopping parameters, i.e., for a given 
Hamiltonian, this yields at most one single solution. Some examples of the wave functions 
for different values of f3 are presented in Fig. ^ 



B. The case # 

In order to obtain the eigenstates described above, we introduced parameters in the 
Hamiltonian ( |1.1[ ) and determined them by requiring that the ansatz ( |3.2D fulfills the tight- 
binding equations ( p.l| ). In Eq. (|1.1|) , we already included the possibility of site-dependent 
on-site energies e^. In the present case, it is natural to choose the on-site energies Ei according 
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to the vertex type of site i. That is, et = with eight parameters /zi, . . . , /is according to 
the eight vertex types of the Penrose tihng. 

Of course, we can perform the same analysis as above for the more general problem 
— it just amounts to replacing the left-hand side of the first three lines of Eqs. (|3.6| ) by 
{E — ^p)Ai with V = 1,4,6, and in the remaining five lines by {E — ^p)A2 with v = 
2, 3, 5, 7, 8, respectively. We do not show the explicit solution of the full problem because it 
is rather lengthy. Although the general solution contains a few free parameters, for a given 
Hamiltonian we still find at most one exact eigenstate. 

In order to compare with Sutherland's resultjll we consider the case without hopping 
along the diagonals of the rhombi, i.e., di = d2 = = = 0. We can express the solutions 
in terms of the three parameters E, f], and 7 := ^42/^1: 





= E - 47 


yW2 


= E 


-2(3- 




= E-2/3-1 -7-1 


/i4 


= E 


— 57 




= E-5(3- 27-1 




= E 


-37 




= E-5(3--f-^ 


1^8 


= E 


-5/3 



(3.16) 

Setting 7 = 1 one recovers Sutherland's solution.0 

Taking into account that we have introduced eight parameters in our the Hamilto- 
nian, it was almost obvious that solutions exist. It would be, however, more interesting to 
introduce additional parameters in the ansatz for the wave function. In this way, one might 
perhaps be able to obtain several eigenstates of a given Hamiltonian and thus come closer 
to the general solution of our problem. This is the subject of the following section. 



IV. GENERALIZED ANSATZ FOR THE EIGENFUNCTIONS 

The Penrose tiling possesses a so-called inflation/ deflation symmetry.il0 In an inflation 
step, the two types of rhombi are dissected into smaller pieces that again constitute a 
rhombic Penrose tiling, but on a smaller scale with all lengths divided by the golden ratio 
r = (1 + 75) /2. The inverse procedure, in which tiles are combined to form larger tiles, is 
known as deflation. 

The idea now is to generalize the ansatz ( |3.2| ) for the wave function by using the vertex 
types and potentials of the deflated tiling in addition to those of the original tiling. Even 
more general, one may consider a sequence of n tilings obtained by successive deflation steps, 
probing the original tiling on larger and larger length scales. In this way, we assign to each 
vertex i of the original tiling a sequence of integers {^kii)}) k = 0,1,..., n, where I'kii) 
specifies the corresponding vertex type in the /c-fold deflated tiling, with k = referring to 
the original tiling. This leads to the following generalized ansatz for the wave function 

n 

#^ = A^M^)} n f^r^"^ (4.1) 

fc=0 

where rrik denotes the double-arrow potential in the /c-fold deflated tiling, and Pk are n + 1 
free parameters. 

It is not completely obvious how to assign the vertex type Ukii) of a site i in the fc-fold 
deflated tiling. Here, we decided to use the concept of the Voronoi cell. We are looking for 
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a Voronoi cell of the deflated tiling that covers the Voronoi cell of our site i in the original 
tiling completely, or at least its largest part. In Fig. ^ we show how the Voronoi cells of the 
original and the two-fold deflated tiling relate to each other. If a cell of the original tiling 
is shared between several larger cells, we assign the vertex to the cell with the maximum 
overlap. However, there are still ambiguities when overlaps of equal area occur. For instance, 
let us concentrate on the case n = 2. In the example shown in Fig. ^ one recognizes that 
the cell corresponding to vertex type 1 (cf. Fig. ^) may be dissected equally between the 
cells corresponding to vertex types 2 or 3 of the deflated tiling. In this case, we cannot 
assign the deflated vertex type unequivocally. Therefore, we demand that the corresponding 
terms in the ansatz ( |4.1| ) are equal. In our example, this yields the equation ^221 = ^321 
for the amplitudes A,j^y^y^ in the ansatz ( [4.1|) , labeled by three digits according to the three 
vertex types. Considering also the flrst deflation step, not shown in Fig. |^, one flnds another 
condition A222 = ^232- 

We now use the ansatz ( [4.1| ) to flnd solutions of the tight-binding equations. Here, we 
restrict ourselves to the case n = 2. In order to set up the equations, we need to consider 
larger patches which may be obtained by two-fold inflation of the 31 second-order vertex 
types of Fig. ^. Each of these patches then leads to a number of equations. Of the 8^ 
possible combinations of indices i/2,i^i,i^0) only 25 occur in the Penrose tiling. Altogether 
we have to deal with a system of 97 equations in 32 variables, namely 25 amplitudes A^^i,-^^^, 
three variables ^2,^1, Po, the four hopping parameters di, ^2, d^, ^4, and the energy E. We 
used Mathematical to solve this system. As above, we flnd three sets of solutions, which 
we express in terms of (32o '■= P2P0 and Pi. They have the following form. 
Solution (1'): The 25 amplitudes, without normalization, are 



^221 


= ^321 


= — /520 (2/^20 ~ 


■ Pi) Po 




^222 


= ^322 


= ^232 = ^532 


= ^732 


— ^832 — P20 Pi Po 


^233 


= ^333 


= P20 Po 






^123 


= ^423 


= Aq23 = ^133 


= ^633 


— P20 


^654 


= y4i74 


= A^sA = —020 


i^Pio - 


-Pi) 


^365 


= ^217 


= ^548 = ^748 


= ^848 


= P20 Pi Po 


^236 


= ^336 


= -{2Pl-Pl 


)P'o 





(4.2) 

and the transfer integrals and the energy are given by the same expressions ( p.9|) as for 
solution (1), where now 

P20P1 ^43^ 



^ Pio ± P i 

In contrast to the amplitudes, the transfer integrals and the energy are hence expressed 
exclusively in terms of P20 and Pi, that is, the hopping parameters and the energy depend 
on P2 and Pq only via the product P2o- 
Solution (2'): Here, the amphtudes read 

^221 = ^321 = P20 Po 

A222 = As22 = A232 = ^532 = ^732 = ^832 = P20 Pi Po 
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^233 


= ^333 


= P20 








^123 


= ^423 


= ^623 = 


^133 


— ^633 


~ P20 


^654 


= Aij4 


= ^484 = 








^365 


= ^217 


= ^548 = 


^748 


= ^848 


= PlPl 


^236 


= ^336 


= P20 Pq 









(4.4) 



and the parameters now follow from the expressions ( p.ll| ) for solution (2) with given by 
Eq. (IJ) and 



:= 4 ; (4.5) 



0i 

thus again they depend on [3q and [32 only via /32o- 
Solution (3'): Finally, 



^221 


= ^321 


= P20 Pi Po 






^222 


= ^322 


= A232 = ^532 = 


^732 


— ^832 — P20 Pi Po 


^233 


= ^333 


= P20 Po 






^123 


= ^423 


= ^623 = ^133 = 


^633 


— P20 


^654 


= Aij4 


= ^484 = P20 Pi 






^365 


= A.217 


= ^548 = ^748 = 


^848 


= P20P1 Po 


^236 


= ^336 


= P'iP'o 







(4.6) 



where again the transfer integrals and the energy follow from the previous expressions (|3.13|) 
for solution (3) by replacing b+ by Eq. ( |4.3| ) and P by Eq. (^3|). 

These solutions comprise those found in the previous section. Indeed, setting P2 = Pi = ^ 
{P20 = Po), we recover the solutions (|3.8| )- (|3.13| ), apart from a common normalization factor 
Pq in the amplitudes. In addition, Eqs. ( |4.2| )-(|4l^) show that the corresponding energy 
eigenvalues are infinitely degenerate. For given values of /?2o and Pi, the Hamiltonian and 
the energy E are fixed, but the eigenf unctions still involve the free parameter j3o- In other 
words, each choice of P2 and Po with the same product yields an eigenstate to the same 
eigenvalue. We note that infinite degeneracies in the spectrum were previously observed in 
tight-binding models on the Penrose tilings. One example is given by the confined degenerate 
states located at the energy E = in the vertex model with di = d2 = ds = (1^ = 0.00 Also 
some of the critical, self-similar eigenstates found in the center model appear to be infinitely 
degenerate.^ 

It is a question whether a larger number of deflation steps, i.e., a larger value of n 
in the ansatz (|4.1j ), leads to further solutions of the tight-binding equations. The larger 
n, the larger is the number of sequences {i'k}o<k<n that occur, and hence the number of 
independent amplitudes. Indeed, for n = 2 we had 25 sequences, for n = 3 and n = 4 
there are 49 and 104, respectively. One might suspect that in the limit n 00, when the 
quantity of sequences tends to infinity, every site is uniquely determined by its sequence, 
and hence one should arrive at the complete solution in the limiting case. However, this is 
not the case, which follows from the fact that — looking at it from the opposite point of 
view — the dissection of a cell under inflation may contain several copies of the same cell 
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type. Therefore, it is doubtful whether larger values of n will lead to new wave functions. 
For n < 4, no solutions beyond ( |4.2D -( ^^ were found. Nevertheless, this does not prove 
that further generalizations might not be more rewarding. 



V. MULTIFRACTAL ANALYSIS 

Already a glimpse at Fig. |^ gives the impression that the wave functions are self-similar. 
Let us therefore investigate this property more thoroughly. To do this, we have to understand 
the distribution of the double-arrow potential m on the tiling. Sutherland^ considered the 
transformation of single and double arrows under two-fold inflation, and proved that the 
value of the double-arrow potential changes at most by 21 under a 2Z-fold inflation. 

For definiteness, let us consider a vertex of type 8 which has double arrows pointing 
outwards in all five directions. In Fig. we show this patch together with its two-fold 
inflation. For the original patch, the values of the double-arrow potential are at the center 
by our choice of normalization, and 1 elsewhere. In the inflated version, the potential takes 
values between and 3, compare Fig. |^. In what follows, we use 2/-fold inflations of this 
particular patch for the multifractal analysis. In this case, the values of the double-arrow 
potential grow linearly with the number of inflation steps. This may be different if one starts 
from other initial patches, for example, starting from vertex type 4 results in a decreasing 
double-arrow potential, corresponding to a different choice of the reference point for the 
potential in the infinite tiling. We note that Refs. ^ and ^ used vertex type 4, together 



with the opposite direction of the arrows, which then also gives an increasing potential. 
Following Refs. ^ and we define a partition function for the 2Z-fold inflated system 

r(^>-5 20 := ^E|£, (5-1) 

where is the norm of the wave function on the finite patch, i.e., N"^ = J2i |</'iP- Here, Si 
denotes the area of the Voronoi cell of vertex i. For a given q, there exists a certain number 
uj{q) such that the partition function ( |5.1| ) is bounded (from above and below) in the limit 
/ — > oo, i.e., it neither vanishes nor diverges. The generalized dimensions 



D, := ^ (5.2) 
g- 1 



completely describe the multifractal properties of the wave function (j). 

In an inflation step, the edge lengths of the rhombi are scaled by a factor r~^. Therefore, 
the area 5*^,^'^ of a Voronoi cell corresponding to vertex type u of the 2/-fold inflated tiling is 
given by 

For simplicity, we restrict our analysis to the ansatz (|3.2|) for the wave function. In fact, 
since only the absolute values of the wave function amplitudes enter in Eq. ( |5.1j ), this also 
applies to the solutions (|3.15|) . Substituting the ansatz into Eq. (|0|) yields 
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T{q,uj-2l) 



m(i) 1 2q 



E 

i/=i 



E 

1/1 |2<j 2/a; 21 

^1 ^ |/5|2^-K(m;20 



(5.4) 



where Viy{m; 21) denotes the number of vertices of type u with potential m multiplied by the 
area of the corresponding Voronoi cell after 21 inflation steps. 

In order to calculate V^(m; 21), we consider the transformation of the Voronoi cells of the 
eight vertex types under a two-fold inflation, compare Figs. ^ and |^. From this, one derives 
recursion relations for the distributions Vy{m\2l) by counting the number of inflated cells 
that are covered by the original cell. For example, as shown in the lower right corner of 
Fig. 1^, the Voronoi cell corresponding to the vertex type 8 with a potential m turns into: (i) 
one cell of type 8 with potential m; (ii) five cells of type 2 with potential m + 1] and (iii) five 
fractional parts (each with an area fraction of (4 — r)/ll ~ 0.216542) of type-6 cells with 
potential m + 1. Conversely, a cell of type 8 in the inflated patch may stem from a vertex of 
type 5, 7, or 8, each of those yielding precisely one complete cell of type 8. Considering all 
vertex types, and computing the fractional areas involved, one arrives at recursion relations 

K(m;2/ + 2) = ^ ^ M^^J F^(m + a; 2/) (5.5) 

cr= — 1 11=1 

with three 8x8 matrices M'^~^\ M^'^\ and M^^\ The quantities we need are certain 
transforms V^(/9; 21) of V^(m; 2Z), defined as 



21 



K(/?;20 := Y.rVu{m-2l) 



(5.6) 



m=0 



see Eq. ( p.4|) . From the recursion relations (|5.5| ), one finds that the transforms Vy{f3\ 21) for 
two successive inflation steps are related by 



K(/?;2/ + 2) = ^M,,^(/5)V;(/5;20 



where the matrix M{(3) reads as follows 



M(/3) 



/ 



21+6T a 
59 

-16+20 



29 
140-6T 



-/3 



52-2T 
59 

95-IOt 



31 
1 








2/3-1 



7+T , 8-2t o-l 
11 ^ 11 I-' 

1 





52-2T g-l 35+IOt ^ 



/3- 
2^/3- 



59 
33-5T 



-20+25T g 153-IOto-I 



77+6T a~l 
31 ^ 





1 

14+2r 
11 



1/3-1 



-12+15T r. 
29 

125-12T 
31 

1 







149-5T o-l 
on P 





5/3-1 






(5.7) 



20-5to_1 20-5t q-1 
11 P 11 P 



(5.8) 
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compare Ref. |T^. It is related to the matrices M^'^^ ( ^.5|) by 

1 

M(/5) = ^ p-^M^^^ 



(5.9) 



cr=-l 



hence the elements of M^'^^ are nothing but the coefficients of of the elements of M{(3). 

The asymptotic behavior (for / — oo) of V^(/3; 21) is governed by the eigenvalue f2(/3) of 
M(/3) with largest modulus 



K(/5;2/) ~ fi'(/5)/.(/5) 



(5.10) 



where /(/3) is the corresponding eigenvector. For positive values of /3, the largest eigenvalue 
(in absolute value) is positive and unique, because the third power of M(/5) is a positive 
matrix. Calculating the norm 



El^.PK(|/3p;2/) 



u=l 



(5.11) 



and substituting the asymptotic behavior of V,y{l3; 21) into ( p^) 



8 1/1 \2q 2luj 

r(g, 21) = y: ,[,::L,... v.m"'-^ 2/) 



^1 (5^°^)"/'Ar29 



r^"l](|/?P'?) 
fi^(|/5P) 



(5.12) 



leads us to the conclusion that the partition function r(g,a;; 21) can be bounded if and only 
if 



uj{q) 



1 



21ogr 



log 



mm 



(5.13) 



In Fig. |10|, we present the fractal exponent Dg ( |5.2|) for several values of (3. For /5 = 1, the 



wave function does not depend on the potential m and takes at most four different values 
according to the translation class of the site. In this case Dg is constant. The smaller 
the faster the wave function decays, leading to a steeper curve Dg as a function of q. 

Concerning the matrix M[j3) ( |5.8| ), we remark that its eigenvalues and eigenvectors are 
connected to the frequencies of the vertex types, which count how often a certain vertex 
type occurs in the Penrose tiling. Indeed, if we set /3 = 1, we obtain a substitution matrix 
for the inflation rules in the Penrose tiling. Therefore, according to the Perron-Frobenius 
theorem, the eigenvector /(I) corresponding to the the eigenvalue with largest modulus f2(l) 
should reproduce the relative frequencies of the vertex types in the tiling. We calculated 
numerically /(I) and found perfect agreement with the known frequencies.SS 

We note that the multifractal analysis can be carried out for the generalized eigenstates 
( [4.2D -( ^^ analogously. However, it becomes more complicated because we have to consider 
the substitution matrix of vertices labeled by the inflated vertex types, which results in a 
25 X 25 matrix. 
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VI. CONCLUSIONS 



We constructed exact non-normalizable eigenfunctions for certain vertex-type tight- 
binding models on the rhombic Penrose tihng. The construction is based on a potential 
m derived from the matching rules of the Penrose tiling that had been introduced in a sim- 
ilar context previously.0ii We consider several generalizations of the ansatz for the wave 
functions, which show that the eigenstates we found are infinitely degenerate. Further 
generalizations can be investigated in a systematic way, and may lead to a wider class of 
accessible wave functions. We hope to report on this, and on the application of this ansatz 
to other quasiperiodic tight-binding models (particularly for the three-dimensional case) in 
the future. 

From the ansatz, it is apparent that the eigenfunctions ( p.8[ ), ( |3.10D and ( |3.12D reflect 
the distribution of the potential m on the lattice. The multifractal analysis of the eigen- 
states therefore reduces to the analysis of the distribution of the potential which was already 
considered by Sutherland.^ It shows that the wave functions are critical, i.e., neither ex- 
tended nor exponentially localized, as typically expected in two-dimensional quasiperiodic 
tight-binding models. 



The present work is a generalization of the ideas of Refs. and and we recover the 
solutions found by Sutherland^ as a special case. In Ref. the authors considered a center 
model on the Penrose tiling, where they also found infinitely degenerate critical eigenstates. 
It is interesting to note that all exactly known eigenstates in such models, including the 
confined states,cl Si appear at energies with infinite degeneracy. At present, we do not 
know whether there is a deeper reason for this observation. 
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FIGURES 




1 1 

FIG. 1. The two types of rhombi in the Penrose tihng and the assignment of hopping integrals 
di, d2, ds, and ^4 to their diagonals. The hopping integral along the edges is chosen as 1. 
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FIG. 3. The double-arrow potential for a patch of the Penrose tiling. 
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FIG. 5. Second-order vertex types corresponding to a central vertex of type 1. Here, the 
encircled numbers denote the vertex types, not the potential. 
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FIG. 6. Wave functions (|3.8| ) for four different values /3 = 0.1, 0.2, 0.6, and 0.9. On a finite 
patch of = 16 757 vertices obtained by sevenfold inflation of a vertex of type 4, the wave 
function has been normalized to J2i \ = 1- The radii of the circles encode in units of the 
edge length they have been chosen as = for 7^|(/>p < 10"^^ R = 1.81og(10W|^n/log(10'^) for 
10-2 < J^\(j)\2 < 102 and R = 1.8 for AA|(/.|2 > 102. 
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FIG. 7. Relation of the Voronoi cells of the Penrose tiling (thin lines) and their two-fold deflation 
(thick lines). The eight patches are obtained by a two- fold inflation of the eight vertex types (see 
Fig. The numbers denote the change of the double-arrow potential with respect to the central 
Voronoi cell. 
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FIG. 8. Voronoi cells of a patch of the Penrose tiling (thin lines) and of its two-fold deflation 
(thick lines). Shaded cells corresponding to vertex type 1 cannot be uniquely assigned to a cell of 
the deflated tiling. 
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FIG. 9. A vertex of type 8 (grey) together with its two- fold inflation (black). 
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